Gap solitons in Bose-Einstein condensates in linear and nonlinear optical lattices 
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Properties of localized states on array of BEC confined to a potential, representing superposition 
of linear and nonlinear optical lattices are investigated. For a shallow lattice case the coupled mode 
system has been derived. The modulational instability of nonlinear plane waves is analyzed. We 
revealed new types of gap solitons and studied their stability. For the first time a moving soliton 
solution has been found. Analytical predictions are confirmed by numerical simulations of the 
Gross-Pitaevskii equation with jointly acting linear and nonlinear periodic potentials. 



I. INTRODUCTION 

The Bose-Einstein condensate in a linear periodic po- 
tential attracts a great attention for last years. Many fas- 
cinating phenomena like Josephson oscillations, macro- 
scopic quantum tunnelling and localization, gap solitons 
etc have been predicted and observed in the experiments 
P, 0| • The description of these phenomena is based on 
the Gross-Pitaevskii equation governing the condensate 
wave function and having terms corresponding to the ex- 
ternal linear (varying in space and time) potential and a 
mean field nonlinearity(taking into account many body 
effects). The strength of the mean field nonlinearity is 
proportional to the atomic scattering length. We can 
imagine the BEC system, when the strength of two body 
interaction is varying in space [1, 0, H, @|- For a peri- 
odic variation of scattering length it leads to appearance 
of a nonlinear periodic potential. The ground state and 
dynamics of localized states of BEC under action of the 
nonlinear periodic potential have been studied recently 
in papers [3, [1, [1] • 

Periodic potential in BEC[l(|[ll| can be generated by 
optical methods. Two types of optical lattices are consid- 
ered: a linear periodic potential induced by the standing 
laser field [l[ and a nonlinear optical lattice produced 
by two counter propagating laser beams with parame- 
ters near the optically induced Feshbach resonance (FR) 
0, Periodic variation of the laser field intensity in 
space by proper choice of the resonance detuning can 
lead to the spatial dependence of the scattering length 
[Tl| . Then the nonlinear term in the GP equation be- 
comes periodically modulated in space, leading to the 
generation of nonlinear optical lattice. 

In real physical situations exact resonance condition is 
not attained and the potential represents the mixture of 
linear and nonlinear optical lattices. This type of periodic 
potential also can be produced by two pairs of counter 
propagating laser beams, when one pair produces linear 
optical lattice, and second pair - a nonlinear one. An- 
other interesting system where the superposition of linear 
and nonlinear periodic potentials can be realized is the 
array of fermion-boson mixtures [12j. Gap solitons(GS) 
in a photorefractive crystal, considered recently in [l3l |. 
also originate from the jointly acting linear and nonlinear 



periodic potentials. But the last system is different from 
ones considered here, since its nonlinearity has saturable 
nature. 

Two types of optical lattices can be distinguished: a 
shallow and a deep. The first type is realized for V < Er, 
and the second one for V > 5Er, where V is the strength 
of the periodic potential and En = h 2 k 2 /2m is the recoil 
energy, k is the laser beam wavenumber. 

In the case of a shallow optical lattice (that will be 
considered in this work) we can apply a coupled mode 
approach. The coupled-mode theory describes general 
properties of gap solitons, taking into account interac- 
tion of counter propagating waves. Analytical solutions 
have been found for the case of the Kerr nonlinearity 
in [l4| . In our case the system of coupled- mode theory 
has more complicated character and is analogous to the 
studied in the nonlinear optics for a deep Bragg grating 
case [HI and a nonlinear layered structures (Tq|. In the 
case of vanishing nonlinear lattice a periodic potential 
transforms into the standard linear one 
study the modulational instability(MI) of nonlinear plane 
waves which is important mechanism for the soliton and 
solitonic pattern generation. The gap soliton solutions 
of this modified coupled mode system are obtained and 
their stability is analyzed. 

The emphasis is also given to the travelling gap soli- 
ton (GS) solutions. We obtain solutions for moving GS 
and study collisions between solitons having an inelas- 
tic character. In optics they have been experimentally 
observed in fibers with Bragg grating[20j. It will be in- 
teresting to observe the travelling GS in BEC with linear 
and nonlinear optical lattices. 

The case of deep optical lattice can be investigated in 
the framework of the tight-binding approximation. The 
GP equation then transforms into the modified discrete 
nonlinear Schrodin ger equation (NLSE) with a nonlocal 
nonlinearity [2H I22I |23| . Analysis of MI and the discrete 
breathers for this model requires separate investigation. 

The paper is organized as follows. In Section II we 
describe the model and derive the coupled mode system. 
In Section III modulational instability of cw solutions is 
investigated. The gap soliton solutions and their stabil- 
ity is discussed in Section IV. The moving gap soliton 
solutions and the solitons collision are studied in Section 
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V. Main results of the paper are summarized in Section 
VI. 



II. THE MODEL. COUPLED-MODE 
EQUATIONS 

We consider here the dynamics of BEC under jointly 
acting linear and nonlinear optical lattices. First type of 
potential is induced by a periodic interference pattern 
from the counter propagating laser beams [Ij. Second 
type of potential is produced by the optically induced 
Feshbach resonance [Tljj . According to this technique, 
periodic variation of laser field intensity in the standing 
wave I(x) = Iocos 2 (kx) produces periodic variation of 
the atomic scattering length a s , i.e. 



&s0 



l 8 + aV 



(1) 



where 8 is the detuning. At large detunings from FR 
a s = a S Q + a s i cos 2 (kx). 

In the recent experiment with the 87 Rb condensate 
the scattering length has been optically manipulated [Ic| . 
The number of of atoms was 10 6 and laser power was ap- 
proximately bOOW/cm 2 . By choice of the laser power and 
detuning of the laser beam around the photo-association 
resonance the scattering length was changed over one or- 
der of magnitude, from 10ao to 190ao (ao is the Bohr 
radius) . 

The quasi ID GP equation for the BEC wavefunction 
under action of superposition of a linear and nonlinear 
optical lattice is 0, [|| : 



iHljj T = - — IpxX 

2m 



V Q cos 2 (kx)t(; + (g {0) 



ID 



+ 5 «cos 2 (fc a; ))M 2 V, (2) 



Collecting terms at e ±lx we obtain the system of coupled 
mode equations 



iA t + 2iA x + l0 {\A\ 2 A + 2\B\ 2 A) 
±B + ^(\B\ 2 B + 2\A\ 2 B + A 2 B*) = 0, 
iB t -2iB x + l0 (\B\ 2 B + 2\A\ 2 B) 



-A 



: (\A\ 2 A- 



2\B\ 2 A- 



B 2 A*) = 0. 



(5) 
(6) 



In comparison with the standard coupled-mode sys- 
tem, this system includes terms corresponding to the 
nonlinear coupling ~ k. 

Note that at small amplitudes of waves A, B ~ 
exp(ikx — iu>t) the dispersion law is co 2 = e 2 /4 + fc 2 , and 
this system has the forbidden gap for — e/2 < u < e/2. 
The system is unchanged if we perform transformations 
7o -» -7o, k -» k, (A, B, u) -> (A\ —B*,—uj). 

It means that results obtained for the first excited 
band u> > e/2 with repulsive interactions can be valid 
for the zero order allowed band lu < —e/2 with attrac- 
tive interactions [18[. 

The Eqs. ©, © read as 



iA(B) t 



SH 



SA*(B*) 

The Hamiltonian H has the following form 
H = J dx[~(A*A x - AA* X - B*B X + BB* 



e -{BA* + AB*) + ^(4\A\ 2 \B\ 2 



\A\ 4 



\B\ 



[(\A\ 2 + \B\ 2 )(BA*+B*A)}. (7) 



The number of atoms is conserved 



TV = 



dx[\A\< 



\B\ 



(8) 



where ^j*^ 1 ^ = 2Ha S Q t iU>±, and a s is the atomic scattering 
length. Introducing dimcnsionlcss variables 



To 



Vo 
' 2Er 

Er 



21,2 



E R = 



K 2 k 

2m 



T 

x = Xk, t = — 
T 



'2|a s0 |wj_ 



Er 



ijje 



a s i 

dsO 



we can rewrite the equation in the form 

iu t + u xx + (70 + Kcos(2o;))|u| 2 u — ecos(2a;)it = 0. (3) 

70 > corresponds to the attractive and 70 < - to 
the repulsive condensates respectively. We will consider 
the shallow lattice case, when e -C 1. In this case one 
can derive the system of coupled mode equations. Let us 
represent the field u(x, t) in the form 



u(x, t) = (A(x, t)e ix + B(x, t)e- ix )e 



(4) 



III. MODULATION AL INSTABILITY 

We start with consideration of the modulational insta- 
bility(MI) phenomena in this model. The MI is a general 
phenomenon originating from the interplay between non- 
linearity and dispersion and it consists in the instability 
of nonlinear plane wave solutions against long wavelength 
modulations. In periodic media the periodicity leads to 
modification of the dispersion law. MI is known to be an 
important mechanism for the generation of solitons and 
solitonic trains in nonlinear optics and Bose-Einstein con- 
densates [H, [25], [2(| . The MI in photonic crystals with a 
deep Bragg grating recently has been investigated numer- 
ically in Ref. [27j . The cw solutions of the coupled- mode 
system are sought using the parametrization suggested 
in [H 



A = 



_ e i(Qx-nt) q _ 



j{Qx-nt) 



(9) 



3 



where a 2 = \A\" 



370 2 e 1 + / 
n = "4 — 



-B| 2 and we have 

2 rea 2 1 + 6/ 2 + f 



Q 



7o« 2 l-/ 2 

4 1 + / 2 



e 1 



4 
/ 2 



/(I + P) 
rea 2 1 - f 2 



f 
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, (10) 
(11) 



When re = the result for a standard coupled-mode sys- 
tem is reproduced. The physical mean of the parameter 
/: / < corresponds to the top of a band gap of the 
linear dispersion curve, and / > to the bottom of a 
band gap . 

To study MI we look for the solution of the form 



where 



N = (e-2K) 2 -2fG{e-2K), 



M a = [2K + 3fG)(e-2K). 

f = 1 and / = — 1 correspond to the bottom (negative 
effective mass) and the top (positive effective mass) of 
the band gap, respectively, The coefficients of Eq. fl6]) 
are real and this equation can be solved analytically 



N 



iV 2 + 16(N + 2M )q 2 



(17) 



A = 



D 



+ 5A(x,t) e ^ Qx - Qt \ 



-^L= + SB(x,t)j e'W-nO, 



(12) 



where 5 A and SB are perturbations of the cw solu- 
tion. Substituting these expressions into Eqs.© and 
using a linear approximation, we get the system for 
5A(x,t),5B(x,t) 



iSA t + 2iSA x + -fSA 



-SB 



K 



SB* + 25B 



-SA" 



G[SA + SA* + 2f(SB + SB*)} 
/(W 2 ) 2/ 

1 + / 2 M+ l + / 2 
i5B t - 2iSB x + ^jSB - ^-SA 

G[f 2 (5B* + SB) + 2f(5A + SA*)} - 



0. 



K 



SA* + 25 A ~ 



I-/ 2 
/(1 + / 2 ) 



SB 



2/ 



1 + / 2 



;6B* 



= 0.(13) 



Here G = j a 2 /(l + / 2 ), K = rea 2 /2. The solution of 
Eqs. (fT3|) can be sought in the form 



8A(B) = C(D) cos{qx - ut) + iE(F) sm(qx - Lit). (14) 

Substituting these expressions into the system (|13[) , we 
get the characteristic determinant for it. The full expres- 
sion is cumbersome and only some numerical simulations 
has been performed for the photonic crystal case[27|. 

Here we perform the analytical consideration, for the 
case / = ±1. 

For Eqs. (fT3|) we obtain the following characteristic 
equation 



U>2 



No + 8q 2 - yV 2 + 16(iV + 2M )<z 2 



(18) 



Complex roots of Eqs. (| 1 T[> . ([18]) correspond to the 
MI of cw solutions. Analysis of these equations shows, 
existence of complex roots is possible in two cases: 

q 2 > -Bo, R < 
and 

q 2 < M , Mo > 
where 



Rq 



N 2 



16(JV + 2M Q ) 



Let us analyze these MI conditions depending on the non- 
linear coupling coefficient, re. There is a threshold value 
re = e/a 2 , when iV = 0, M = 0, R = and the MI 
condition q 2 > — i? is replaced by q 2 > M . 
We consider two cases: 

a) / = 1 The root cj 2 is complex at 

re > re , q 2 > — i?o; « < «0i <? 2 < M . 

b) / = — 1 The root 101 is complex at 

re > re + 70, q 2 > -Ro; 



ko + 7o < « < «0j 



re < re g 2 > — i?o- 

As to the root u>i it is complex for all values of re at 
q 2 > —Ro except two points re = reo + 70 and re = kq. 

When / = 1 and re = the instability of cw exists in 
the interval 



uj a - [(e - 2^) 2 - 2fG(e - 2K) + 8q 2 } lu 2 - 

8q 2 (2K + 3fG){e-2K) + I6q 4 = 0. (15) 



It should be reminded that here |/| = 1. 
Let us rewrite this equation in the form 

lu A - (8q 2 + N )lu 2 - 8M a q 2 + Wq 4 = 



(16) 



^eG<q< J^eG, 



that is the same as the result obtained for the standard 
coupled mode system (26[. 

The dependence of MI gain versus the wavenumber of 
modulations of plane wave is presented in Fig. [TJ The 
dashed line corresponds to 7 = 1, a = 1 and solid line to 
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FIG. 1: MI gain versus the wavenumber of modulations q. 
Parameters are k = 1, e = 0.2, / = — 1. Solid line corre- 
sponds to the case 70 = 0, a = 0.35, and dashed line to 
70 = 1, a = 1 . 



FIG. 3: Suppressed modulational instability. Parameters are 
70 = 0, k = 1, e = 0.2, a = 0.8, / = -1 q = 0.1 . 




FIG. 2: Modulational instability development. Parameters 
are 70 = 1, ft = 1, e = 0.2, a = 1, / = — 1. a) g = 0.1, b) 
q = 0.5 



7 = 0, cv = 0.35. Other parameters of the system are the 
same in both cases: e = 0.2, k = 1, / = — 1. 

Results of numerical simulations of MI when 
70 = 1,k = 1, e = 0.2, a = 1,/ = — 1 is presented 
in Fig. O The evolution of amplitudes for q = 0.1 and 
g = 0.5 are shown. Initial value of the wave function 
amplitude at t = is 1.44. By the time t = 20 the 
amplitude value grows up to 5.83. 

A curious behaviour of the modulational instability has 
been observed when numerical simulating MI at 7 = 0. 
(Other parameters of the system were k = 1, e = 0.2 and 
initial cw amplitude a — 0.8). The result of simulation 



FIG. 4: MI gain versus the initial wavefunction amplitude 
60^- Parameters are k = 1, e = 0.2, / = — 1. Solid line 
corresponds to the case 70 = 0, a = 0.8 and dashed line to 
70 = 1, a = 1 . 



is shown in Fig. [3] 

In this case one can observe suppressed MI evolution 
with no growth in the wave function amplitude. Such a 
development of MI can be explained from the dependence 
of MI gain versus the cw amplitude a given in Fig. 3] 
As seen, in the case of 7 = MI gain exists in some 
limited range 0.42 < a < 0.82. Taking into account that 
a = max \^\ and that outside of this limited range the 
solution is stable, we can conclude that the growth of the 
wave function amplitude must be suppressed. 

In general dependence of MI gain on the system param- 
eters is very complicated that is demonstrated in Figs. [S] 
- [7J Results of numerical simulations confirm predictions 
of the theory. 
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FIG. 5: MI gain versus the wavenumber of modulations q 
and the initial wavefunction amplitude a. Parameters are 
k = 1, e — 0.2, f — — 1, 70 = — 1. (The same is for / = 
1, 70 = 1.) 



FIG. 8: Phase-plane portrait U = U(9, H g ) at 70 = 1, e = 0.2, 
K = 1 and = —0.04. Solid line is for H — 0, dashed line for 
= -0.04, dotted line for H = 0.00085. 



IV. GAP SOLITON IN LINEAR-NONLINEAR 
LATTICE 




FIG. 6: MI gain versus the wavenumber of modulations q 
and the initial wavefunction amplitude a. Parameters are 
k = 1.1, e = 0.2, f = —1, 70 = 1. (The same is for / = 
1, 70 = -1.) 




FIG. 7: MI gain versus the wavenumber of modulations q 
and the initial wavefunction amplitude a. Parameters are 
k — 1, e = 0.2, f — —1, 70 = 1. (The same is for / = 
1, 70 = -1.) 



To find gap soliton solution let us look for the partic- 
ular solution of the form A, B — > (a(x), b(x)) exp(— iuot). 

The symmetry of the system ([5]) allows to consider the 
case a(x) = b* (a:)[28j. Then we get the equation for a(x) 



uja + 2ia x + 37o|a| a — -a* + 



— K\a\ 2 a* 
2 1 1 



~a d = 0. 



(19) 



The solution will be sought in the form 

12 



a 



Substituting this expression into Ea. ([Tn)) we get the fol- 
lowing system of equations 

Q x = -u - 3j Q U+ I cos 9 - 2kU cos 9, (20) 

U x = |Z7 sin 61- nU 2 sm9. (21) 

The system has the Hamiltonian H g and can be written 
in the form 



dHg 
dU 



Ux ~ 89 ' 



where 



H g = -luU - t^JoU 2 + |[/cos0- kU 2 cost 



(22) 



(23) 



We will consider solutions that rapidly decay in the 
space, i.e. U — > 0, |x| — * 00. In this case H g = 0, that 
corresponds to separatrix in the phase-plane portrait (see 
Pig. H. 

Then for 9 we have the following equation 



■ cos 9, 



(24) 



6 



with U determined as 

U = 

Integrating Eq. 

cos 9 
where 



2uj — e cos ( 



(25) 



370 + 2k cos 9' 
provided that | tan((9/2)| < 7 we find 
1 - 7 2 tanh 2 (/?x) 



1 + 7 2 tanh (f3x) 



(26) 



7 = 



e - 2to 



Ve 2 - 4w 2 



e + 2w ' 4 
From Eq. (|25[) we obtain expression for the amplitude 



U 

U(x) 



e-2uj 



(370 + 2k) cosh 2 (/3x) + (37o - 2k) 7 2 sinh 2 (/3x) 
Finally the gap solution takes the form 



u(x, t) — 2\fij cos i^x 



9(x) 



-i(u+l)t 



with 



- e/2 < uj < e/2. 
The number of atoms in the soliton is 

/oc 
|u(x,t)| 2 dx 
-OC 

When 370 > 2k 

8 

AT = 



Udx. 



and when 370 < 2k 



: arctan 



(370 - 2k)7 
V97o 2 - 4« 2 ' 



v/4« 2 - 9 7o 2 
In the case 



: In 



3 7o + 2k + 7 ^4k 2 - 9 7o 2 



3 7o + 2k - 7v /4k 2 - 9 7 2 



370 > 2k 



(27) 
(28) 

(29) 
(30) 

(31) 

(32) 
(33) 



solution P?| is stable for all values of ui within condition 
(|2^|) . i.e. the soliton exists everywhere in the gap. 

In the case 370 < 2k for stable solutions we have an- 
other condition for u) 



3e7o e 
— — < co < -. 
4k 2 



(34) 



Chemical potential oj versus norm N is depicted in 
FigEl According to the Vakhitov-Kolokolov criterion [29] 
the soliton is stable provided that dN/du < 0. 

Evolution of the stable solution for 7o = 1, e = 0.2 
is given in Fig. 1101 The initial wave packet envelope is 
taken in the form of Eq. (|2~T|) . 



0.1 



-0.1 




FIG. 9: The chemical potential ui versus norm TV. The curve 
(a) corresponds to condition (|33|) and the curve (b) to (|34[) . 




FIG. 10: Stable solution for 70 = 1, e = 0.2, k = 1 and 
uj = 0.05 at different times : (a) t = and (b) t = 100. The 
initial wave packet envelope is taken in the form of Eq. (|27[) . 



An interesting case we have when 7o = 0. It can be 
realized with the Feshbach resonance method. Then the 
soliton exists in the interval < lu < e/2, i.e. in the upper 
half of the gap. The number of atoms for uj = A <C e 
near the center of the gap can be estimated as 



N 



7T A 



When lu — e/2 — A, i.e. the parameter near the band 
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FIG. 11: Stable solution for 70 = 0, e = 0.2, k = 1 and 
w = 0.05 at different times : (a) t = and (b) t = 100. The 
initial wave packet envelope is taken in the form Eg. (|27|l . 



edge, the number of atoms in gap soliton tends to zero 
like 



TV 



4 A 

7T V £ 

Numerical simulation of this case is shown in Fig. 1111 

We performed numerical simulations of the GP equa- 
tion with the linear-nonlinear optical potentials. In Figs. 
[TTil [TT1 evolution of the initial state in Eq. © is given in 
the form of Eq. ([27j) . 

When I tan(0/2)| > 7 the expressions for the amplitude 
U and cos 9 take the forms 



COS* 



tanh 2 (/?x) - 7 2 
tanh 2 (/3a;) + j 2 



(35) 



U(x) 



(370 + 2k) smh z (Px) + (370 - 2k) 7 2 cosh 2 (fix) 

In this case the solution is evidently unstable. Indeed, 
in this case dN /dui > and in accordance with the 
Vakhitov-Kolokolov criterion 2!l such a solution is un- 
stable. PDE simulation of the evolution of an unstable 
solution for 70 = 1, e = 0.2 and k — 2 is given in Fig. [T2l 
The initial wave packet envelope is taken in the form of 
Eq.®. 



V. TRAVELLING GAP-SOLITON SOLUTION. 
SOLITONS COLLISION 

Travelling solitons can be obtained proceeding from 
the anzatz of work[l5j. The solution of Eqs.([5][H]) should 
be sought in the form of the sum of travelling waves with 
different amplitudes and phases 



A(x,t) 
B(x,t) - 



i(M£)-f)t) 



(37) 
(38) 




300 



FIG. 12: Unstable solution for 70 = 1, e = 0.2, k = 2 and 
u = -0.08 at different times : (a) t = 0, (b) t = 40 and (c) 
t = 140. The initial wave packet envelope is taken in the form 
Eq.(T 



where 



£ = x — vt. 



Substituting these expressions into Eqs.([S][H]) we have 



F 6 = -Fsin(6 B -6 A )- HF 2 sm(0 B -6 A ), (39) 



370F - 2kF co$(6 b - 9a), 



•a), -J 



noF 



4 + V 
kFcos(9b — 9 a 



(40) 
(41) 



where F = Fa = F B 



7o 



(4 + v 2 ) 



A 3 7o, A 



Ae, k = A k, 
1 



(42) 



The first two equations of this system coincide with 
Eqs. (|20M2ip with renormalized parameters. Thus the 
region of the solitons stability are determined again by 
conditions Eqs. (f2"9"]) , (f3~4"|) with renormalized parameters 
as in Eq. (|42]) . 

The gap is given by 



e v 2 e v 2 

- TTVl r < ^ < -\ 1 -. 

2 V 4 2 V 4 



(43) 



When the velocity tends to the limiting value v — > 2, the 
gap region shrinks to zero. Also the gap does not exist 
for the velocity v > 2. 
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400 



FIG. 13: The time dependence of the travelling soliton speed 
for jo = 1, e = 0.2, k = 1 and uj = 0.05 when the initial speed 
v = 0.2. 




FIG. 14: Collision of two moving solitons for 70 = 1, e = 0.2, 
k, = 1 and cj = 0.05. 



elastic. A small radiation at the solitons tails seems to 
be caused by the nonintegrability of the system. 



VI. CONCLUSION 

In conclusion, we have investigated the dynamics of 
BEC under joint action of the linear and nonlinear op- 
tical lattices. These lattices can be created by counter 
propagating laser beams via optically induced Feshbach 
resonances. Also such structures can be generated in ar- 
rays of the Fermi-Bose mixtures. For the case of shallow 
lattices we derive the coupled mode system of equations 
and have investigated the modulational instability of cw 
matter waves in such periodic media. We find the insta- 
bility regions of cw nonlinear matter waves. The knowl- 
edge of MI regions can be useful for generation of gap 
solitons in such structures. New types of gap soliton so- 
lutions are found and the stability regions are analyzed. 
In particular, the gap soliton can exist in the upper half 
of the forbidden gap, when the background atomic scat- 
tering length is equal to zero(7o = 0). We find the trav- 
elling gap soliton solution in BEC and study collision of 
gap solitons by numerical simulations. It is shown that 
the collisions are nearly elastic, with the emitting of small 
radiation, due to the nonintegrability of the couple-mode 
system. 
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The stability condition is 

3(4 + v 2 )ej ~ : , ■ 

; — < <^ < ^\ 1 ^ ~r- ( 44 ) 

4k 2 V 4 v ; 

Numerical simulations with the initial wave packet 
taken in the form 

u{x, t) = {A(Qe ix + B(£)e lx )e- l{u+1)t (45) 

show that obtained solution of Eq. ([3]) has the form of an 
automodel soliton. 

In the course of evolution of the initial wave-packet 
the soliton parameters, the soliton speed, amplitude and 
width, approach their stationary values. Evolution of the 
soliton speed is shown in Fig. [T3J 

To check the nature of obtained moving solutions in 
Fig. Q3] the result of simulation of two solitary waves 
collision taken in the form of Ea.(|45p is shown. As is 
evident from this figure, collision of two solitons is nearly 
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